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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  inversion  approaches  that  enable  the  estimation  of 
refractivity  profdes  and  the  associated  uncertainty.  Further  to  develop  methods  for  mapping  the 
refractivity  parameters  and  their  associated  uncertainty  into  propagation. 

OBJECTIVES 

The  development  of  new  parametric  approaches  for  the  inversion  of  radar  clutter  data  to  estimate 
atmospheric  refractivity  including  an  assessment  of  the  quality  of  the  parameter  estimates. 

APPROACH 

Our  inversion  approach  has  mainly  been  based  on  SAGA  [Gerstoft  97]  and  focused  on  estimation  of 
the  parameters  corresponding  to  the  field  that  gives  the  best  fit  to  the  data.  We  have  concentrated  on 
demonstrating  the  feasibility  of  RFC  using  an  efficient  1 1 -parameter  description  of  the  environment. 
The  quality  of  the  inversion  was  addressed  by  comparing  the  field  using  the  estimated  parameters  to  a 
measured  field  [Gerstoft  03a,  Gerstoft  03b,  Gerstoft  03c].  Little  has  been  done  to  indicate  the  quality 
of  the  solution  for  each  parameter,  either  with  the  variance  of  parameter-estimate  or  preferably  the 
complete  a  posteriori  distribution.  We  have  already  done  much  work  on  this  in  an  ocean  acoustic 
context  [Gerstoft  98],  but  this  has  not  been  explored  in  our  RFC  processing  to  date.  This  will  entail 
developing  likelihood  formulations  and  importance  sampling  algorithms.  This  inversion  approach  will 
show  the  information  content  in  the  data,  the  importance  of  each  parameter,  and  the  quality  of  the 
inversions. 

WORK  COMPLETED 
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An  algorithm  for  estimation  of  propagation  loss  and  its  statistical  properties  based  on  radar  sea  clutter 
has  been  developed  [Gerstoft  03d]. 

RESULTS 

The  simulated  data  are  generated  based  on  the  helicopter  measured  range-dependent  refractivity 
profiles  (Run  7)  for  the  Wallops  98  experiment  [Gerstoft03a],  see  Fig  1.  A  range  interval  from  10-100 
km  is  used.  A  simple  trilinear  model  (Fig  2)  is  used  for  the  inversion  of  the  refractivity  profile  as 
outlined  in  the  Appendix  of  [Gerstoft03a].  We  then  search  for  tri-linear  refractivity  parameters  at  0 
and  100  km  range.  To  obtain  refractivity  profiles  at  other  ranges  the  parameters  are  interpolated 
linearly.  The  first  3  parameters  were  given  a  uniform  distribution  but  the  slope  was  given  a  non- 
uniform  distribution  as  indicated  in  Fig  3  (left  bottom).  This  is  because  a  negative  mixed  layer  slope  is 
only  likely  for  shallow  ducts  (low  base  height). 

First,  90,000  environmental  models  are  selected  from  the  prior  distributions,  Fig  3  left.  Using  these 
environmental  models,  90,000  vectors  of  propagation  loss  versus  range  are  precomputed.  The 
posteriori  distribution  (right  column  Fig.  3)  is  computed  as  a  product  of  the  prior  distributions  with  the 
likelihood  distribution.  The  likelihood  distribution  shows  how  well  the  environmental  models  describe 
the  data.  We  note  that  the  overall  behavior  of  the  posteriori  distributions  seems  reasonable. 


2 


O  1 00  200  300  400 


0  1 00  200  300  400 


Baseheight  (m) 


Baseheight  (m) 


2 


1 

0  50  100 

M— deficit  (M) 


2 


1 

O  50  100 

Thickness  (m) 


Slope  (M/m) 


Figure  3  Prior  (left)  and  posteriori  (right)  distributions  of  the  parameter  estimates.  The  abscissa  “1” 
and  “ 2  ”  refers  to  the  values  at  emitter  and  receiver,  respectively. 

The  overall  objective  is  being  able  to  estimate  posteriori  statistics  of  propagation  loss.  An  example  of 
this  is  shown  in  Fig  4.  We  compute  the  average  priori  propagation  loss  based  on  an  even  weighting  of 
the  propagation  loss  from  each  generated  refractivity  model.  The  posteriori  probability  distribution  of 
the  propagation  loss  is  based  on  weighting  the  propagation  loss  from  each  refractivity  model  with  the 
posteriori  probability. 

The  average  prior  and  posteriori  propagation  losses  are  shown  in  Fig.  4b  and  c.  It  is  seen  that  the 
average  posteriori  propagation  loss  identifies  a  ducting  environment  as  observed  in  the  data,  Fig.  4a, 
but  the  prior  does  not. 

The  probability  distribution  of  the  propagation  loss  then  is  computed  at  all  ranges  and  depths.  Both 
prior  (Fig.  5  left)  and  posteriori  (Fig.  5  right)  propagation  loss  are  shown  at  10  and  100  m  height  as  a 
function  of  range.  Clearly,  the  posteriori  plot  shows  how  the  data  has  improved  our  estimation  of 
propagation  loss. 
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Figure  4:  The  propagation-loss  field  from  the  true  environment  (from  Fig.  la)  and  the  average 
propagation-loss  field  based  on  (b)  prior  information,  and  (c)  posteriori  information  (bottom). 
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Figure  5  Prior  (left)  and  posteriori  (right)  propagation  loss  probability  distributions  at  10  m  (top) 

and  100  m  (bottom). 


IMPACT/APPLICATIONS 

Knowledge  of  refractivity  profiles  is  important  for  radar  performance  prediction.  A  very  promising 

approach  is  to  determining  refractivity  based  on  the  radar  clutter  return. 

RELATED  PROJECTS 

Refactivity  Data  Fusion  and  Assimilation  (Ted  Rogers,  SPAWAR):  This  project  is  concerned  with  near 

real-time  techniques  for  inferring  refractivity  parameters  from  radar  sea  clutter. 
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